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and  the  corresponding  electron  density  profiles  generated  from  the 
SSIES/SM  data  for  the  8  January  1988  pass.  The  x-axis  labels  in  (a)  are 

for  the  satellite  location,  those  in  (b)  are  for  the  3(X)-km  field-line 
"footprint."  The  x-axis  labels  are  GLAT  for  geographic  latitude,  GLON  for 
geographic  longitude,  ALT  for  satellite  altitude,  APXLAT  for  modified 
apex  latitude,  APXLON  for  apex  longitude,  and  APXLT  for  apex  local 
time. 

2  Data  from  Figure  1(a)  replotted  to  the  "window"  sampled  by  the  E1SCAT  8 

radar.  Note  the  reversal  of  north  and  south  from  Figure  l(u)  to  this  figure. 

3  Data  from  the  8  January  1988  campaign  day:  (a)  detrended  intensity  data  16 

from  the  A10-AFSAT  VHF  radio  link,  (b)  detrended  phase  data  from  this 

link,  and  (c)  total  ion  density  from  the  SSIES/SM  sensor.  The  plots  are 
aligned  to  cover  roughly  the  same  spatial  region. 

4  Comparison  of  S4  (a)  and  (b)  scintillation  indices  calculated  from  the  18 

AIO-AFSAT  intensity  and  phase  data  (solid  lines)  to  the  corresponding 
indices  calculated  from  the  SSIES  C*L  estimates  for  two  irregularity 
models  (17:2  wings  plotted  as  squares,  20:1  rods  as  diamonds)  and  those 
calculated  from  the  WBMOD  model  (crosses).  The  two  points  in  (a) 
labeled  "retune"  included  retune  discontinuities,  and  the  point  labeled 
"focus"  contained  a  strong  focus.  The  S4  values  for  these  points  are  unre¬ 
liable. 

5  Total  ion  density  measured  by  the  SSIES/SM  sensor  on  DMSP  F8  near  20 

Kwujalein  on  24  August  1988.  The  y-uxis  is  ion  density,  ranging  from  0.0 

to  2.0x10s;  and  the  x-axis  labels  are  GMT  for  Greenwich  Mean  Time 
(HH:MM:SS),  APXLAT  for  modified  apex  latitude,  and  APXl.uN  for  apex 
longitude. 
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PREFACE 


This  report  summarizes  the  work  completed  during  an  investigation  into  the  feasibility  of 
using  in-situ  observations  of  the  ionosphere  from  the  DMSP  SS1ES  sensors  to  calculate  parame¬ 
ters  that  characterize  ionospheric  scintillation  effects.  This  work  is  part  of  a  larger  effort  with  an 
overall  objective  of  providing  the  USAF  Air  Weather  Service  with  the  capability  of  observing 
ionospheric  scintillations,  and  the  plasma  density  irregularities  that  cause  the  scintillations,  in 
near  real-time  and  updating  models  of  ionospheric  scintillation  with  these  observations. 


VI 


1.  Introduction 


Many  modern  military  systems  used  for  communications,  command  and  control,  naviga¬ 
tion,  and  surveillance  depend  on  reliable  and  relatively  noise-free  transmission  of  radiowave 
signals  through  the  earth's  ionosphere.  Small-scale  irregularities  in  the  ionospheric  density  can 
cause  severe  distortion,  known  as  radiowave  scintillation,  of  both  the  amplitude  and  phase  of 
these  signals,  A  basic  tool  used  in  estimating  these  effects  on  systems  is  a  computer  program, 
WBMOD,  based  on  a  single-scatter  phase-screen  propagation  model  and  a  number  of  empirical 
models  of  the  global  morphology  of  ionospheric  density  irregularities! 1  *al.  An  inherent  weakness 
of  WBMOD  is  that  the  irregularity  models  provide  median  estimates  for  parameters  with  large 
dynamic  ranges,  which  can  lead  to  large  under-  and  over-estimation  of  the  effects  of  the  iono¬ 
spheric  irregularities  on  a  system, 

One  solution  to  this  problem,  at  least  for  near  real-time  estimates,  is  to  update  the 
WBMOD  irregularity  models  with  observations  of  the  various  parameters  modeled.  One 
proposed  source  for  these  observations  is  from  the  in-situ  plasma  density  monitor  to  be  flown  on 
the  Defense  Meteorology  Satellite  Program  (DMSP)  satellites.  Previous  studies!3'  using  in-situ 
measurements  from  the  DE-2  satellite  have  found  that  there  is  potential  for  using  the  data  in  this 
fashion.  This  study  is  designed  to  assess  the  applicability  of  this  data  set  to  real-time  updates  of 
the  WBMOD  models.  There  are  two  primary  objectives: 

(1)  Develop  and  refine  techniques  for  generating  estimates  of  parameters  that  char¬ 
acterize  ionospheric  scintillation  from  in-situ  observations  of  the  ionospheric  plasma  from  the 
DMSP  SSIES  sensors. 

(2)  Determine  if  the  parameters  calculated  from  the  SS1ES  data  can  be  used  to 
compute  the  scintillation  effects  on  a  transionospheric  radiowave  signal. 

This  report  describes  one  final  aspect  of  the  analysis  of  data  from  a  campaign  conducted 
near  Trotnso  in  January  1 9KH  and  summarizes  the  findings  and  results  of  the  entire  project. 
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2.  Background 

The  propagation  model  used  in  the  WBMOD  program  (based  on  weak-scatter  phase- 
screen  theory!'-2!)  characterizes  the  ionospheric  electron  density  irregularities  that  cause  scin¬ 
tillation  via  eight  independent  parameters!4!; 

(1)  a:  The  irregularity  axial  ratio  along  the  direction  of  the  ambient  geomagnetic  field. 

(2)  b:  The  irregularity  axial  ratio  perpendicular  to  the  direction  of  the  ambient  geomagnetic 
field. 

(3)  5;  The  angle  between  sheet-like  irregularity  structures  and  geomagnetic  L  shells. 

(4)  hp:  The  height  of  the  equivalent  phase  screen  above  the  earth's  surface, 

(5)  iq:  The  in-situ  irregularity  drift  velocity. 

(6)  a0:  The  outer  scale  of  the  irregularity  spectrum . 

(7)  q:  The  slope  of  a  power-law  distribution  that  describes  the  one-dimensional  power- 
density  spectrum  (PDS)  of  the  irregularities. 

(8)  C*L:  The  height-integrated  strength  parameter. 

The  first  three  parameters  (a,  b,  and  8)  and  the  direction  of  the  ambient  geomagnetic  field  specify 
the  propagation  geometry,  while  the  last  three  (a0,  q,  and  C*L)  specify  the  spectral  characteris¬ 
tics  of  the  irregularities. 

It  may  be  possible  to  obtain  estimates  for  the  values  of  three  of  those  parameters  from  the 
DMSP  SSIES  sensors:  (from  the  SSIES  Ion  Drift  Meter  (DM)),  and  q  and  QL  (from  the 
SSIES  Ion  Scintillation  Meter  (SM)).  In  this  study,  we  will  focus  on  the  estimation  of  QL  from 
this  data  set  and  consider  q  and  only  in  terms  of  the  effects  of  uncertainties  in  these  parame¬ 
ters  on  the  estimates  of  C*L.  Of  the  eight  parameters,  C^L  varies  the  most  as  a  function  of  loca¬ 
tion  and  time,  and  has  the  most  profound  effect  on  the  accuracy  of  estimates  of  scintillation 
levels  made  by  the  WBMOD  model. 

In  the  phase-screen  propagation  theory  used  in  WBMODPI,  the  C*L  parameter  is  actually 
the  product  of  two  parameters:  Q ,  the  three-dimensional  spectral  "strength"  of  the  electron 
density  irregularities  at  a  scale  size  of  1  km*  (related  to  the  structure  constant  used  in  classical 
turbulence  theory);  and  L,  the  thickness  of  the  irregularity  layer.  The  models  in  WBMOD  were 
obtained  from  analysis  of  phase  scintillation  data  from  the  WIDEBAND  and  HiLat  satellites, 
which  will  provide  estimates  of  the  height-integrated  value  of  C^L  rather  than  independent 


*  The  cited  reference  develops  the  theory  in  terms  of  an  earlier  definition  of  the  strength  param¬ 
eter,  Cs,  which  is  defined  at  a  scale  size  of  2 n  meters,  It  is  related  to  C*  according  to  the  equa¬ 
tion  Cs  =  (2it/l  000)4+2  ck  . 
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measures  of  Cj.  and  L,  Because  of  this,  the  model  was  developed  for  Q.L  rather  than  for  and 
L  separately. 


The  calculation  of  an  estimate  of  the  C*L  parameter  from  topside  in-situ  ion  density 
observations  requires  two  operations.  First,  an  estimate  of  C*  at  the  satellite  altitude  is  made 
from  a  finite-length  time  series  of  density  measurements.  Second,  the  estimate  of  C*  is 
converted  to  an  estimate  of  C*L  in  some  fashion  which  will  account  for  both  the  thickness  of  the 
irregularity  layer  and  the  variation  of  C* ,  or  <ANC2>,  within  the  layer. 

The  data  set  from  which  the  estimates  of  these  parameters  are  to  be  obtained  will  be 
collected  by  three  instruments  in  the  DMSP  SSIES  (Special  Sensor  for  Ions,  Electrons,  and 
Scintillation)  sensor  package.  This  data  set  will  contain  the  following  in-situ  observations: 

(1)  High  time-resolution  (24  samples/sec)  measurements  of  the  ion  density  and 
measurements  of  the  ion  density  irregularity  PDS  at  high  fluctuation  frequencies  from  the  Ion 
Scintillation  Meter  (SM)l5', 

(2)  Measurements  of  the  horizontal  und  vertical  cross-track  ion  drift  velocities  from  the 
Ion  Drift  Meter  (DM)I*I. 

(3)  Measurements  of  the  ion  and  electron  temperatures,  the  densities  of  0+  and  the 
dominant  light  ion  (H+  or  He+),  and  the  horizontal  rum  ion  drift  velocity  from  the  ion  Retarding 
Potential  Analyzer  (RPA)i6!, 

The  basic  data  of  this  set  are  the  high  time-resolution  density  data  from  the  SM  which  will  be 
used  to  generate  estimates  of  the  irregularity  PDS.  The  drift-velocity  measurements  from  the 
DM  and  RPA  will  be  used  in  calculating  an  estimate  of  C*  from  parameters  obtained  from  the 
PDS,  and  the  other  measurements  from  the  RPA  will  be  used  in  calculating  C\L  from  C* . 

In  the  first  year  of  this  project,  techniques  for  calculating  estimates  of  C*L  from  the 
SSIES  data  set  were  developed,  and  parametric  studies  were  conducted  to  determine  the  uncer¬ 
tainties  in  the  final  C*L  estimates  due  to  uncertainties  in  the  parameters  and  procedures  used  to 
calculate  the  estimates.  Since  the  first  SSIES  sensor  puckage  was  not  flown  until  mid- 1987, 
these  studies  were  conducted  using  simulated  SM  density  data  sets  and  phase  scintillation  data 
from  the  Wideband  satellite,  The  results  of  these  studies  were  reported  in  Scientific  Report  1  for 
this  project!7!  (herein  referred  to  as  Repoit  1)  and  will  be  summarized  in  Section  4.1  of  this 
report. 


The  second  phase  of  the  project,  which  focuses  on  how  well  these  techniques  work  using 
data  from  the  DMSP  F8  and  F9  spacecraft,  was  begun  during  the  second  year  of  the  project.  A 
coordinated,  multi-sensor  observation  campaign  was  conducted  during  January  1988  in  the 
vicinity  of  Tromso,  Norway,  in  order  to  collect  data  for  this  study.  The  GL  Airborne  Iono¬ 
spheric  Observatory  (AlO)  aircraft  flew  repeated  north-south  legs  along  the  magnetic  meridian  at 
times  when  the  DMSP  F8  satellite  was  passing  overhead.  Data  collected  on  the  AIO  included 
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intensity  and  phase  scintillation  observations  on  a  VHF  link  from  an  AFSAT  satellite,  auroral 
images  from  an  all-sky  photometer,  and  ionograms  from  a  digisonde.  Ionospheric  soundings  by 
the  EISCAT  incoherent- scatter  radar  located  in  Norway  also  were  made  nperiodically  throughout 
the  observation  period,  Data  were  collected  for  DMSP  passes  on  8, 9, 15,  16,  17,  and  18  January 
1988,  The  procedures  used  to  process  these  data  and  results  from  the  analysis  were  reported  in 
Scientific  Reports  2  and  3  for  this  project!8'9!  (herein  referred  to  as  Repons  2  and  3)  and  will  be 
summarized  in  Section  4.2, 

During  the  third  year,  a  complete  analysis  was  made  SSIES  data  collected  from  both 
DMSP  F8  and  F9  during  a  multi-sensor  campaign  conducted  by  the  Defense  Nuclear  Agency 
(DNA)  near  Kwajalein  Island  were  analyzed  to  determine  the  usefulness  of  these  data  for  scin¬ 
tillation  characterization  in  the  equatorial  ionosphere,  The  results  of  these  studies  were 
presented  in  Report  3  and  urc  summarized  in  Section  4.3  of  this  report. 


A  note  on  coordinate  systems;  The  geomagnetic  latitude,  longitude,  and  local  time  coor¬ 
dinate  system  used  throughout  this  report  is  modified  apex  coordinates! l°l,  a  coordinate  system 
derived  from  apex  coordinates  proposed  by  VanZundt  et  a/1"!.  This  system,  which  is  used  in  the 
WBMOD  model,  was  chosen  because  it  is  very  similar  to  both  invariant  and  corrected 
geomagnetic  coordinates  at  high  latitudes,  and  modified  apex  latitude  is  very  similar  to  dip  lati¬ 
tude  in  the  equatorial  region,  becoming  identical  to  dip  latitude  at  the  dip  equator. 
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3.  Analysis  of  SSIES/AIO/EISCAT  Campaign  Data 

In  Report  3,  the  results  from  a  detailed  analysis  of  the  data  collected  on  8  January  1988 
were  presented.  One  aspect  of  the  data  analysis  which  was  not  presented  in  detail  was  a  discus¬ 
sion  of  the  electron  density  profiles  used  to  generate  the  effective  layer  thickness  parameter, 
Le  f  f ,  which  is  used  in  turn  to  calculate  an  estimate  of  C*L  from  the  value  derived  from  the 
in-situ  Scintillation  Meier  (SM)  ion  density  measurements.  In  this  section  we  will  describe  the 
way  in  which  these  profiles  were  generated  and  present  examples  of  the  profiles  used  in  the  8 
January  analysis. 

The  model  used  to  generate  the  profiles  is  a  two-component  diffusive-equilibrium  model 
developed  several  years  ago  for  building  profiles  from  SSIES  obseivations  (see  Appendix). 
Inputs  to  the  model  include  the  critical  frequency  of  the  F2  peak  and  the  height  of  the  peak  (^2 
and  hmF2),  the  semi-thickness  of  the  topside  just  above  the  F2  peak,  the  electron  density  (or  total 
ion  density)  at  some  height  in  the  topside,  and  either  the  height  at  which  the  0+  and  H+  number 
densities  are  equal  (i.e.,  the  O+/H+  transition  height)  or  the  ratio  of  0+  to  H+  at  some  height  in 
the  topside,  In  addition,  the  model  has  several  "free"  parameters  which  can  be  adjusted  to 
change  the  shape  of  the  profile.  Any  of  the  input  parameters  can  be  adjusted  within  the  model  to 
fit  the  profile  to  the  other  inputs.  For  example,  in  the  present  study  the  f0F2  and  hmF2  values 
were  adjusted  to  fit  the  profile  to  the  total  ion  density  observed  by  the  SM  sensor. 

The  procedures  to  build  the  profile  used  to  calculate  Lef  f  were  as  follows: 

(1)  Calculate  the  average  total  ion  density  for  the  data  set  used  to  calculate  C* 
from  the  density-trend  time  series  removed  from  the  SM  density  measure¬ 
ments  in  the  detrending  process. 

(2)  Trace  down  the  field  line  passing  through  the  center  of  the  data  sample  to  an 
altitude  of  300  km. 

(3)  Calculate  estimates  of  f0F2  and  hmF2  using  the  CCIR  coefficients  for  f0F2 
and  M3000  and  values  for  the  0+/H+  transition  height  and  the  topside  semi¬ 
thickness  at  this  location. 

(4)  Iteratively  adjust  the  f„F2  and  hmF2  parameters  to  fit  the  profile  to  the  aver¬ 
age  density  calculated  in  step  (1). 

This  profile  and  a  model  for  the  height  variation  of  the  irregularity  strength  as  a  function  of  the 
background  plasma  density  are  then  used  to  calculate  Leff  > 

Figure  1  shows  the  SM  density  record  from  the  8  January  1988  pass  near  Tromso  (upper 
plot)  and  a  contour  plot  of  the  corresponding  electron  density  profiles  which  were  generated  to 
calculate  Leff:.  The  locations  and  geomagnetic  local  time  (APXLT)  along  the  x-axis  of  the  in- 
situ  density  plot  are  calculated  at  the  satellite  altitude,  and  those  along  the  x-axis  of  the  profile 
plot  are  calculated  at  the  3()()km  "footprint"  of  the  field  line  passing  through  the  satellite  location. 
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Figure  1.  Total  ion  density  data  from  the  SSIES  Scintillation  Meter  (SSIES/SM)  (a)  and  the 
corresponding  electron  density  profiles  generated  from  the  SSIES/SM  data  for  the  8 
January  1988  pass.  The  x-axis  labels  in  (a)  are  for  the  satellite  location,  those  in  (b)  are 
for  the  300- km  field-line  "footprint."  The  x-axis  labels  are  GLAT  for  geographic  latitude, 
GLON  for  geographic  longitude,  ALT  for  satellite  altitude,  APXLAT  for  modified  apex 
latitude,  APXLON  for  apex  longitude,  and  APXLT  for  apex  local  time. 
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The  slight  tilt  to  the  right  in  the  profile  plot  is  due  to  mapping  the  profiles  constructed  along  field 
lines.  (Note:  The  profile  plot  shows  no  ionization  in  the  E  region  because  the  model  used  docs 
not  include  an  auroral  precipitation-produced  E  region,  and  there  is  no  solar  component  as  the  E 
region  is  dark  throughout  the  pass.  This  does  not  affect  the  results,  however,  as  the  model  used 
to  calculate  C*L  assumes  that  all  contributions  to  the  integrated  irregularity  strength  come  from 
the  F2  region.] 

Our  initial  hope  was  to  use  profiles  from  the  EISCAT  radar  as  part  of  the  study,  but  there 
have  been  difficulties  in  putting  the  radar  data  into  useful  form.  One  figure  showing  a  contour 
plot  of  electron  density  from  a  roughly  north-south  meridional  scan  from  the  radar  with  a  center 
time  of  18:34:40  GMT  (66800  seconds  past  midnight)  has  been  received  for  comparison  with  the 
profiles  used  in  the  Lef  f  calculation.  Unfortunately,  this  figure  is  in  color  and  could  not  be 
legibly  reproduced  in  black  and  white  for  this  document. 

We  have,  however,  replotted  a  subset  of  the  model  profiles  in  the  same  format  as  the 
color  figure  showing  what  the  radar  would  have  seen  if  it  had  sampled  our  model  ionosphere 
(Figure  2),  This  figure  has  the  same  contour  levels  as  in  Figure  1,  but  the  x-axis  has  been 
reversed  (north  is  now  to  the  right)  and  is  plotted  in  the  corrected  geomagnetic  (CG)  latitude  of 
the  3(X)  km  "footprint"  location,  The  irregular  shape  of  the  sampled  region  corresponds  to  that 
shown  in  the  color  plot  provided  of  the  EISCAT  data.  The  model  parameters  used  to  build  these 
profiles  and  the  resulting  Lef  f  values  are  listed  in  Table  1. 

Quantitatively,  the  radar  data  and  model  "prediction"  agree  fairly  well.  Both  show  an  F2- 
region  enhancement  at  the  northern  edge  of  the  sampled  region,  a  lesser  enhancement  at  the 
southern  edge,  a  "tongue"  of  enhanced  density  reaching  from  the  northern  enhancement  to  the 
center  of  the  sampled  region,  and  a  region  of  low  plasma  in  the  topside  between  CG  latitude 
65.5°  and  68°.  There  arc  some  discrepancies  in  details.  The  altitudes  of  the  peak  in  the  model 
profiles  average  about  30-50  km  above  the  estimated  peaks  in  the  radar  datu,  and  the  enhance¬ 
ment  at  the  northern  end  of  the  region  is  further  from  Trornso  in  the  radar  data  than  is  shown  in 
the  model  profiles.  Unfortunately,  no  quantitative  comparison  of  the  density  levels  was  possible 
as  the  color  plot  of  the  radar  data  was  of  relative  rather  than  absolute  density. 
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Figure  2,  Datu  from  Figure  1(a)  replotted  to  the  "window"  sampled  by  the  EISCAT  radar.  Note 
the  reversal  of  north  and  south  from  Figure  1(a)  to  this  figure. 
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Tabic  1 .  Parameters  for  the  electron  density  profiles  plotted  in  Figure  2. 
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4.  Summary  of  Project  Results 

4.1  Study  of  Methods  for  Calculation  of  C*L 

The  goal  of  this  study  was  to  develop  a  method  for  calculating  an  estimate  of  C/L  from 
data  collected  by  the  SS1ES  plasma  probes  and  to  make  a  determination  of  the  uncertainties  in 
these  estimates. 

4.1.1  Procedures  for  Calculating  C*L 

According  to  phase-screen  theoryl2!,  if  the  power-density  spectrum  (PDS)  of  the  plasma 
density  irregularities  are  assumed  to  follow  a  single-regime  power-law  model,  an  estimate  for  the 
spectral  strength  parameter,  C*,  can  be  calculated  from  in-situ  measurements  of  the  plusmu 
density  irregularities  from 


/lO3^1 

Ck  »  5.0x10s  q  ^ - j  T,  [1] 

where  q  and  Tj_  are  the  slope  and  intercept  of  a  log-linear  fit  to  the  PDS  calculated  from  the 
plasma  density  data,  and  vp  is  an  effective  velocity  of  the  satellite  with  respect  to  the  irregulari¬ 
ties.  This  last  parameter  is  defined  in  a  coordinate  system  defined  by  the  geometry  of  the  irreg¬ 
ularities  und  is  given  by 

vp  -  vT  C  _v  [2) 

where  v  is  the  vector  velocity  of  the  satellite  with  respect  to  the  irregularities  (which  includes  the 
satellite's  orbital  velocity  and  the  in-situ  velocity  of  the  irregularities),  and  C  is  a  transformation 
matrix  derived  from  a  generalized  model  of  the  irregularity  geometry I>21. 

The  procedure  developed  for  calculating  parameters  q  and  Tj  from  the  SSIES  Scintilla¬ 
tion  Meter  ion  density  data  is  as  follows; 

(1)  Detrend  the  ion  density  data  by  removing  a  background  trend.  This  trend  can 
be  calculated  from  either  a  low-pass  filter  applied  to  the  entire  data  set  or  a 
quadratic  least-squares  fit  to  each  data  sample  (256  or  512  points). 

(J.)  Apply  a  time-domain  window  to  the  detrended  data  sample.  The  window  of 
choice  is  a  30%  split-bell  cosine  taper. 

(3)  Use  an  FFI’  of  the  detrended,  windowed  data  sample  to  calculate  an  estimate 
of  the  PDS  (0>iN). 
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(4)  Smooth  the  PDS  estimate.  The  smoothing  function  used  is  a  five-point,  cen¬ 
tered  function  with  binomial  weights. 


(5)  Use  a  linear  least-squares  fit  to  log(PDS)  versus  log(frequency)  over  the  fre¬ 
quency  range  0.3  to  5.0  Hz. 

The  various  parameters  in  this  method,  such  as  the  severity  of  the  windowing  function  and  the 
frequency  range  for  the  PDS  fit,  were  determined  from  a  series  of  parametric  studies  using 
simulation  data  generated  with  known  spectral  characteristics  (see  Report  1)  and  were  validated 
ugainst  a  sample  of  data  from  the  SM  sensor  on  DMSP  satellite  F8  (see  Report  2). 

The  calculation  of  and  estimate  of  C*L  from  an  observation  of  C*  at  the  satellite  altitude, 
;ts,  is  based  on  the  following  assumptions: 

(1)  The  altitude  at  which  the  Q  measurement  is  made  is  within  the  irregularity 
layer, 

(2)  The  measured  Q.  is  representative  of  conditions  throughout  the  irregularity 
layer, 

(3)  The  geometry  of  the  irregularities,  and  the  slope  of  the  irregularity  PDS,  are 
relatively  constant  throughout  the  height  range  of  the  irregularity  layer. 

(4)  The  variation  of  C*,  or  <AN2>,  with  altitude  throughout  the  irregularity  layer 
is  known,  or  can  be  modeled. 


These  assumptions  lead  the  following  relationship  between  C\L  and  C*: 

C^l,  =  Ck  l<eff 


[31 


where  Lpf  f  is  an  effective  layer  thickness  parameter,  given  by 

<AN2>  dh 


<AN2>s 

where  hb  is  the  height  of  the  bottom  of  the  irregularity  layer  and  <AN2>S  is  the  density  variance 
at  the  satellite  altitude. 

r\ 

For  this  project,  we  have  chosen  a  simple  model  for  the  variation  of  <AN^>  with  altitude 
which  assumes  that  the  layer  extends  from  1/4  of  a  scale  height  below  the  F2-luyer  peak  height 
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to  the  satellite  altitude,  and  that  the  ratio  of  <AN2>  to  the  background  density  is  constant 
throughout  the  layer.  This  changes  Equation  (4)  to 


Leff 


/ho 

N  2  (h)  dh 
"b 


Ne2(hs) 


15] 


where  Ne  (h)  is  the  height  variation  of  the  electron  density.  We  have,  used  a  two-component 
(0+  and  H+)  model  of  the  topside  ionosphere  developed  specifically  for  use  with  the  SS1ES  data 
set  for  Ne  (h)  l|3>.  The  application  of  this  model  to  the  problem  at  hand  was  described  in  Section 
3. 


In  summary,  estimates  of  CtL  were  calculated  from  the  SSIES  SM  sensor  data  using  the 
following  steps: 

(1)  Calculate  q  and  from  the  SM  ion  density  data  records  using  either  256-  or 
512-point  data  samples. 

(2)  Calculate  the  effective  velocity,  vp,  from  the  satellite  orbital  velocity,  obser¬ 
vations  or  model  values  for  the  In-situ  drift  velocity  of  the  irregularities,  and 
models  for  the  irregularity  geometry  taken  from  the  WBMOD  scintillation 
modeH2!. 

(3)  Calculate  Q.  using  Equation  1 1 1. 

(4)  Fit  the  Ne  (h)  model  to  the  average  plasma  density  of  the  data  sample  by 
adjusting  the  F2  peak  density  and  height  parameters  (sec  Section  3). 

(5)  Calculate  l.ief  f  using  the  fitted  model  and  Equation  |5|. 

(6)  Calculate  C*L  using  Equation  1 3). 

4.1.2  Uncertainties  in  C*l, 

The  effects  of  uncertainties  in  various  parameters  used  in  the  calculation  of  C^l.  on  the 
final  C*L  were  investigated  and  quantified  (sec  Sections  3  and  4  in  Report  1),  Three  sources  of 
uncertainty  studied  were  errors  in  the  values  of  q  and  T  } ,  errors  in  vp.  and  errors  in  L,e  f  t , 

The  errors  inherent  in  the  calculation  of  q  and  from  a  plasma  density  sample  were 
estimated  using  u  set  of  192  simulated  data  samples  generated  with  known  spectrul  characteris¬ 
tics,  including  values  for  q  varying  across  the  nominal  expected  runge  of  this  parameter  (1.0  to 
3,0)  (see  Appendix  A  of  Report  1  for  the  details  of  this  simulated  data  set).  The  RMS  error  in 
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estimating  q  and  Tj  from  this  data  set  using  the  procedures  outlined  earlier  were  3.3%  and 
6.7%,  respectively.  The  resulting  uncertainty  in  C*  (and,  ultimately,  in  C*L)  due  to  these  RMS 
"error  bars"  is  a  function  of  the  effective  velocity  (Vp),  ranging  from  6  to  11%  for  values  of  Vp 
of  l(XX)  to  7000  nVs.  The  average  uncertainty  due  to  this  source  should  be  on  the  order  of  8- 
10%. 


The  effective  velocity  is  a  function  of  the  satellite  velocity  (vs),  the  in-situ  drift  velocity 
of  the  irregularities  (vcl),  und  the  irregularity  geometry  which  includes  the  orientation  of  the 
irregularities  und  their  uxial  rutios  (a  und  b),  The  velocity  of  the  satellite  due  to  its  orbital 
motion  is  fairly  well  determined,  as  is  the  orientation  of  the  irregularities  (determined  to  first 
order  by  the  local  geomagnetic  field  direction),  so  the  major  sources  of  uncertainty  in  this 
parameter  are  in  the  in-situ  drift  velocity  and  the  axial  rutios  of  the  irregularities.  In  order  to 
reduce  the  complexity  of  the  problem,  the  parametric  studies  conducted  were  limited  to  the 
geometries  expected  for  a  satellite  in  a  i  ominal  DMSP  orbit. 

The  effects  of  uncertainties  in  these  parameters  on  C^L  are  fairly  complex,  but  cun  be 
summarized  as  follows; 

(1)  Uncertainties  in  axial  ratio  a  had  a  minimal  effect  on  C^L.  A  ±50%  uncer¬ 
tainty  in  a  results  in  less  than  a  2%  uncertainly  in  C)L  at  high  latitudes,  and 
rarely  in  more  than  5%  at  other  latitudes. 

(2)  Uncertainties  in  axial  ratio  b,  which  specifics  how  elongated  the  irregularities 
are  across  the  geomagnetic  field,  had  a  much  larger  effect  at  high  latitudes 
(there  is  no  effect  at  other  latitudes,  as  shcct-like  irregularities  occur  only  in 
the  auroral  region).  A  ±50%  uncertainty  in  b  can  result  In  a  25-100%  uncer¬ 
tainty  in  C(U. 

(3)  If  measurements  of  the  in-situ  drift  velocity  are  available  from  the  SSIES 
Drift  Meter  and  Retarding  Potential  Analyzer,  the  expected  uncertainties  in 
these  observations  will  lead  to  5-10%  uncertainties  in  C*L.  If  these 
measurements  are  not  available,  the  average  uncertainty  in  Q.L  increases  to 
10-20%  and  may  be  as  high  as  30-40%, 

As  a  function  of  location,  the  uncertainties  in  C*l,  due  to  uncertainties  in  Vp  in  the  equutorial 
regie'll  should  be  on  the  order  of  5  to  10%  when  the  in-situ  drift  velocities  are  known  and  10  to 
20%  when  they  are  not;  in  the  auroral  region  they  should  be  5  to  10%  when  vcj  and  b  are  known, 
15  to  30%  when  v(1  is  not  well  known  (i.e.,  not  observed),  and  25  to  >KX)%  when  neither  are 
well  known. 

The  effects  of  uncertainties  in  Lof  f.  on  Q.U  were  studied  by  varying  parameters  of  the 
density  profile  model  around  values  nominal  for  regions  which  will  be  sampled  by  the  DMSP 
satellites.  The  results  of  this  study  were  as  follows: 
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(1)  The  most  crucial  profile  parameter  for  making  accurate  estimates  of  C*L  is 
the  plasma  density  at  the  F2  peak.  An  error  of  20%  in  f0F2  can  result  in 
errors  of  100%  in  C*L, 

(2)  The  effects  of  errors  in  the  density  at  the  satellite  can  be  as  serious  as  those 
due  to  errors  in  f0F2,  but  this  parameter  will  be  measured  routinely  with  an 
accuracy  of  about  ±5%,  which  results  in  a  ±8-10%  uncertainty  in  C*L. 

(3)  The  uncertainty  in  C*L  is  much  more  sensitive  to  errors  in  the  ratio  of  the 
density  at  the  satellite  to  the  peak  density  than  to  absolute  (correlated)  errors 
in  the  two  values.  It  is  important  that  the  two  densities  be  consistent  with  one 
another. 

(4)  The  effects  of  20%  errors  in  the  other  profile  parameters,  most  of  which 
describe  the  shape  of  the  profile,  resulted  in  less  than  10%  errors  in  C*L. 

4.2  Auroral  Region  Studies 

In  order  to  assess  the  limitations  of  using  QL  calculated  from  the  DMSP  SSIES  data  to 
characterize  transionospheric  scintillation  in  the  auroral  zone,  a  datu-collection  campaign  was 
conducted  in  the  vicinity  of  the  E1SCAT  incoherent  radar  facility  at  Tromso,  Norway,  during 
January  1988.  The  primary  objective  of  the  SSIES/AIO/EISCAT  campaign  was  to  make  near- 
coincident  observations  of  the  state  of  the  ionosphere  and  of  a  radio  signal  propagating  through 
the  ionosphere  in  order  to  assess  the  validity  of  using  information  from  the  .SSIES  SM  sensor  to 
characterize  and  quantize  the  effects  of  scintillation  on  a  transionospheric  radio  transmission. 
The  concept  was  to  use  estimates  of  the  scintillation  parameters  C*L  and  q  derived  from  the 
SSIES  SM  density  data  to  calculate  estimates  of  the  level  of  phase  and  amplitude  scintillation  to 
be  expected  on  u  specific  transionospheric  satellite  communications  link,  and  then  to  compare 
these  to  the  scintillation  levels  actually  observed  on  that  link. 

The  following  data  were  collected  during  this  campaign: 

(1)  From  the  DMSP  F8  satellite: 

(a)  In-situ  observations  of  the  total  ion  density  (Nj)  at  a  datu  rate  of  24 
samples/sec  from  the  SSIES  SM  sensor. 

(b)  In-situ  observations  of  the  horizontal  (u^,)  and  vertical  (uv)  cross-track  com¬ 
ponents  of  the  ion  drift  velocity  at  a  data  rate  of  6  samples/sec  from  the 
SSIES  DM  (Drift  Meter)  sensor. 

(c)  In-situ  observations  of  precipitating  electrons  and  ions  at  a  data  rate  of  one 
energy  spectrum  per  second  from  the  SSJ/4  sensor. 
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(2)  From  the  AFGL  Airborne  Ionospheric  Observatory  (AIO): 


(a) 

Phase  and  amplitude  records  of  a  250MHz  radio  transmission  from  an 
AFSAT  satellite  (#7225). 

V 

(b) 

All-sky  photometer  images  ut  three  wavelengths  (4278A,  5577A,  and 
6300A)  at  a  data  rate  of  (at  least)  one  image  per  wavelength  per  minute 
from  the  All-Sky  Imaging  Photometer  (ASIP). 

(c) 

Ionospheric  soundings  at  a  data-rate  of  approximately  one  sounding  every 
two  to  three  minutes  from  the  Digisonde. 

(3)  From  the  EISCAT  incoherent  radar: 

(a)  Electron  density  profiles  along  a  ncar-meridian-plane  scan, 

(b)  Electron  and  ion  temperature  profiles  along  a  near-meridian-plane  scan, 

(c)  Ion  drift  velocity  profiles  along  a  ncar-mcridian-plane  scan. 

(Note:  We  had  hoped  to  process  data  from  the  SSIES  Retarding  Potential  Analyzer  (RPA) 
sensor  to  obtain  observations  of  the  ion  temperature  and  the  ratio  of  the  number  density  of  0+ 
ions  to  H+  (or  He+)  ions,  both  of  which  would  have  been  used  in  the  calculation  of  C*L. 
Unfortunately,  this  instrument  failed  prior  to  the  first  cumpaign  period.) 

One  data  set  was  collected  during  each  day  of  the  campaign  (8-9  and  15-18  January 
1988),  The  DMSP  F8  satellite  is  in  a  nominal  dawn-dusk  orbit  which,  during  the  two  campuign 
collection  intervals,  pussed  the  satellite  near  the  EISCAT  facility  at  around  1830UT.  The 
observing  plan  was  to  have  the  AIO  fly  a  north-south  pattern  configured  so  that  the  ionospheric 
penetration  point  of  the  AIO-AFSAT  radio  link  passed  neur  the  DMSP  orbit  track,  and  for  the 
EISCAT  facility  to  conduct  north-south  meridian  scans  roughly  aligned  with  the  DMSP  track  ut 
the  time  of  the  overpass.  Estimates  of  C^I.  were  calculated  for  each  of  the  six  passes,  and  the 
results  were  catalogued  in  Appendix  A  of  Report  2,  Due  to  limited  resources  and  time,  the  anal¬ 
ysis  focused  on  the  data  collected  on  8  January  1988.  The  detuils  of  this  analysis  cun  be  found  in 
Section  3  of  Report  3. 

Figure  3  summarizes  the  basic  data  used  in  the  analysis.  Plots  (a)  and  (b)  in  this  figure 
are  the  intensity  and  phase  data,  respectively,  from  the  AIO-AFSAT  VHF  beacon  link,  and  Plot 
(c)  is  the  Scintillation  Meter  (SM)  density  data  from  DMSP  satellite  F8.  The  time  intervul  of  the 
SM  data  plot  was  chosen  so  that  the  geomagnetic  latitudes  covered  by  this  plot  is  identical  to  that 
covered  by  the  AIO-AFSAT  plots.  Estimates  for  C*L  and  q  were  generated  from  the  SM  data 
which  were  then  used  to  calculate  the  scintillution  parameters  S4  and  c0  using  the  phase-screen 
propagation  model  in  the  WBMOD  scintillation  codeO"»i  and  two  models  for  the  irregularity 
axial  ratios  (20:1  rods  and  17:2  wings),  The  same  parameters  were  calculated  (in  the  time 
domain)  from  the  AIO-AFSAT  intensity  and  phase  data,  respectively. 
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Figure  3.  Data  from  the  8  January  1988  campaign  day;  (a)  detrended  intensity  data  from  the 
AIO-AFSAT  VHF  radio  link,  (b)  detrended  phase  data  from  this  link,  und  (c)  total  ion 
density  from  the  SSfES/SM  sensor,  The  plots  are  aligned  to  cover  roughly  the  sume 
spatial  region. 


Figure  4  shows  the  results  of  these  calculations,  The  solid  lines  in  these  figures  are  the 
indices  calculated  from  the  beacon  data  in  Figure  3,  and  the  diamonds  and  squares  are  values 
calculated  from  the  SM  C^L  and  q  values,  The  crosses  indicate  values  generated  by  the 
WBMOD  models  for  C^L  and  q,  The  agreement  between  the  S4  values  generated  from  the 
DMSP  data  and  that  calculated  from  the  beacon  data  is  quite  good,  particularly  in  light  of  the 
many  assumptions  which  went  into  the  C*L  calculation,  The  agreement  between  the  0$  values  is 
almost  as  bad  as  the  agreement  between  the  S4  values  was  good.  The  validity  of  the  phase  data 
has  been  called  into  question,  however,  by  the  presence  of  unexplained  high  levels  of  low- 
frequency  power  in  the  phase  data  which  does  not  appear  to  be  correlated  with  anything 
geophysical.  Attempts  were  made  to  remove  this  low-frequency  power  by  detrending  and  by 
comparing  the  phase  and  in-sitit  power-density  spectra,  but  no  clear,  supportable  results  were 
obtained. 

The  final  conclusion  of  this  study  was  that  the  excellent  agreement  between  the  S4  values 
predicted  by  the  SM  density  data  and  those  observed  on  the  AIO-AFSAT  radio  link  indicate  that 
there  is  potential  for  using  the  SM  data  in  the  auroral  region  for  characterizing  transionospheric 
scintillation  levels,  This  assessment  should  be  confirmed,  however,  by  a  similar  analysis  of  data 
collected  during  u  second  campaign  in  the  Tromso  area  during  January  1990.  In  this  second 
campaign,  AFSAT  beacon  data  were  collected  both  on  the  AIO  aircraft  and  on  a  ground  receiver 
near  Tromsol14!. 

4.3  Equatorial  Region  Studies 

During  the  month  of  August,  1988,  the  Defense  Nuclear  Agency  (DNA)  conducted  an 
extensive  multi-sensor  duta-collection  campaign  neat  Kwajalein  Islund  to  assess  the  effects  of 
ionospheric  scintillation  on  radar  signul  propagation  (Propagation  Effects  Assessment  - 
Kwajalein,  or  PEAK).  This  campaign  was  used  as  a  target  of  opportunity  to  conduct  an  assess¬ 
ment  of  the  potential  for  using  the  DMSP  SSIES  datu  for  characterizing  scintillation  in  the 
equatorial  region,  The  purpose  of  the  DNA  PEAK  campaign  was  to  collect  data  from  a  number 
of  ionospheric  probes,  both  remote  and  in-situ ,  which  could  be  used  to  characterize  plasma 
density  irregularities  in  the  equatorial  ionosphere  in  order  to  assess  the  effect  of  these  irregulari¬ 
ties  on  transionospheric  radar  propagation,  The  campaign  ran  during  the  period  3-31  August 
1988,  with  most  of  the  ground-bused  instruments  located  on  islands  in  the  Kwajnlein  Atoll 
(9°24'N,  167°28'E).  Of  the  various  data  sets  collected,  the  following  have  been  made  uvailable 
for  this  study: 

a,  SSIES  data  from  both  F8  and  F9  DM.SP  satellites  for  the  entire  month  of  August. 

b,  Intensity  scintillation  data  (S4)  from  a  UHF  link  from  Kwajalein  Island  to  a  FLTSAT 
satellite  located  over  the  equator  at  172°  E  longitude  for  the  period  3-31  August. 
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Figure  4.  Comparison  of  S4  (a)  and  a*  (b)  scintillation  indices  calculated  from  the  AIO-AFSAT 
intensity  and  phase  data  (solid  lines)  to  the  corresponding  indices  calculated  from  the 
SSIES  C*L  estimates  for  two  irregularity  models  (17:2  wings  plotted  us  squares,  20:1 
rods  as  diamonds)  and  those  calculated  from  the  WBMOD  model  (crosses).  The  two 
points  in  (a)  labeled  "retune"  included  retune  discontinuities,  and  the  point  labeled 
"focus"  contained  a  strong  focus.  The  S4  values  for  these  points  are  unreliable. 
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c,  Intensity  scintillation  data  (S4)  from  UHF  and  L-band  links  to  the  HiLat  and  Polar 
BEAR  satellites  obtained  at  Kwajalein  Island  with  the  DNA  ROVER  receiver  for  the  period  1-29 
August. 

In  this  study,  attempts  were  made  to  correlate  observations  of  plume  structures  in  the  SSIES  DM 
density  data  (such  as  that  shown  in  Figure  5)  with  the  FLTSAT  and  HiLat/Polar  BEAR  scin¬ 
tillation  observations,  and  to  corroborate  a  threshold  effect  observed  in  DMSP  Retarding  Poten¬ 
tial  Analyzer  data  in  an  earlier  studyl,5>  (this  will  be  denoted  the  RPA  Study).  Details  of  the 
study  can  be  found  in  Section  4  of  Report  3. 

The  results  of  this  study  were  consistent  with  and,  for  the  most  part,  confirmed  those 
reported  in  the  RPA  Study,  The  main  results  were  as  follows: 

(1)  Due  to  the  inclination  of  the  DMSP  orbit,  the  orientation  of  the  equatorial 
plume  structures  with  the  local  geomagnetic  meridian,  and,  possibly,  the 
altitude  of  the  DMSP  orbit,  there  is  a  significant  probability  that  the  satellite 
may  pass  through  an  equatorial  region  populated  by  plume  structures  and  not 
intersect  one.  There  were  several  cases  in  the  study  period  where  S4  on  the 
FLTSAT  link  was  saturated  for  most  of  the  evening  yet  there  were  no  plume- 
structure  signatures  in  the  DMSP  data. 

(2)  If  u  plume-structure  signature  is  found  in  the  DMSP  data,  radio  links  in  the 
same  longitude  sector  will  experience  periods  of  severe  scintillation. 

(3)  There  may  be  a  threshold  effect  in  the  background  density  observations 
which  might  be  used  to  infer  the  presence  of  plume  structures  even  though 
none  were  encountered  directly  by  the  satellite.  Similar  effects  were  seen  in 
both  studies,  but  the  threshold  levels  found  differed  by  as  much  as  a  factor  of 
2,  and  there  may  be  an  unacceptably  high  false-alarm  rate. 


The  conclusion  of  this  study  was  that  the  DMSP  data  would  be  of  limited  use  for  making 
direct  observations  of  plume  structures  in  the  equatorial  night  region,  but  that  the  data  might  be 
used  in  indirect  ways  to  infer  the  probability  of  the  presence  of  these  structures.  A  more  detailed 
study  using  a  multi-year  database  of  DMSP  density  data  and  any  available  equatorial  scintilla¬ 
tion  data  should  be  conducted  to  develop  and  test  such  indirect  methods,  If  there  is  a  density- 
threshold  effect,  the  results  of  an  expanded  study  might  also  provide  a  basis  for  using  the  output 
of  improved  models  of  the  background  plasma  density  as  an  input  to  the  equatorial  model  in 
the  WBMOD  scintillation  model. 
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DMBP  F9  -  24  Aug  19B8  -  East  Pa»» 
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Hgure  5,  Total  ion  density  measured  by  the  SSIES/SM  sensor  on  DMSP  F8  near  Kwujalein  on 
24  August  1988.  The  y-axis  is  ion  density,  ranging  from  0.0  to  2.0x10s;  and  the  x-axis 
labels  are  GMT  for  Greenwich  Mean  Time  (HH:MM:SS),  APXLAT  for  modified  apex 
latitude,  and  APXLON  for  apex  longitude. 


5.  Conclusion 


The  primary  objective  of  this  project  was  to  assess  the  potential  for  using  in-s it u  obser¬ 
vations  of  the  ionospheric  plasma  available  from  the  DMSP  SSIES  sensor  package  to  character¬ 
ize  the  level  of  scintillation  effects  on  transionospheric  radio  paths.  Methods  were  developed  for 
using  the  data  in  this  fashion,  and  estimates  were  made  of  the  accuracy  of  the  parameters  gener¬ 
ated  using  these  methods.  Data  from  a  multi-instrument  data-collection  campaign  conducted 
near  the  EiSCAT  incoherent  iadar  facility  in  the  early-evening  auroral  region  were  analyzed  to 
test  these  methods  and  assess  the  potential  for  their  use.  Data  from  a  second  campaign  in  the 
equatorial  region  were  analyzed  to  assess  the  usefulness  of  the  SSIES  data  in  determining  the 
presence  of  depletion-plume  structures  in  the  post-sunset/pre-midnight  equatorial  ionosphere. 

The  methods  developed  for  calculating  irregularity/scintillation  parameters  from  the 
SSIES  observations  require  information  not  directly  available  from  the  SSIES  data  set,  in  partic¬ 
ular  the  shape  of  the  plasma  density  irregularities  being  sampled  (the  irregularity  axial  ratios) 
and  the  height  variation  of  the  irregularity  strength.  If  these  data  are  to  be  used  operationally,  the 
missing  information  must  be  provided  either  from  direct  observations  or  from  models.  As  it  is 
highly  unlikely  that  observations  of  either  the  irregularity  shapes  or  the  height  variation  of  the 
irregularity  strength  will  be  available  for  operational  use,  models  must  be  used  to  provide  this 
information. 

A  model  for  the  shape  of  the  irregularities  is  available  from  the  WBMOD  scintillation 
model,  and  while  it  provides  reasonable  estimates  for  the  axial  ratio  parameter  a  (specifying  the 
elongation  of  irregularities  along  the  geomagnetic  field),  there  is  some  concern  about  the  esti¬ 
mates  provided  for  axial  ratio  parameter  b  (specifying  the  presence  and  severity  of  elongation 
across  the  geomagnetic  field  direction)!  Iftl.  Parametric  studies!7!  showed  that  the  irregularity- 
parameter  estimates  calculated  from  the  SSIES  data  ure  very  sensitive  to  changes  in  b,  with 
errors  as  high  as  a  motor  of  two  in  cases  where  this  parameter  was  set  to  specify  rod-like  irregu¬ 
larities  (b  =  1)  when  in  fact  cross-field  elonguted  sheet-like  irregularities  were  present  (b  >  1). 
Further  work  should  be  (lone  to  improve  these  models. 

A  simple  model  was  developed  for  the  height  variation  of  the  irregularity  strength 
(AN  (h) ).  This  model  was  based  on  the  assumption  that  the  relative  irregularity  strength  (AN/N) 
is  invariant  with  height  and  uses  a  two-component  diffusive-equilibrium  model  of  the  electron 
density  variation  in  the  topside  ionosphere!  W  to  provide  the  height  variation  of  N.  An  effective 
layer-thickness  parameter  (Lef  f)  is  calculated  from  the  AN  (h)  model  and  is  used  to  settle  the 
insitu  measurement  of  the  irregularity  strength  (C*)  to  an  integrated  measure  (C*L)  which  char¬ 
acterizes  the  entire  layer.  This  scaling  process  is  very  sensitive  to  the  ratio  of  the  plasma  density 
at  the  satellite  to  that  at  the  F2  layer  peak,  and  less  so  to  details  of  the  height  variation  in  the 
topside.  As  there  will  be  a  fair  measure  of  the  density  at  the  satellite  from  the  SSIES  data,  it  is 
imperative  that  an  accurate  measure  (observation  or  model)  of  the  F'2  peak  density  consistent 
with  the  density  at  the  satellite  be  available. 
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The  results  of  the  high-latitude  assessment  campaign  were  positive,  but  they  were  not 
conclusive  due  to  problems  with  the  phase  scintillation  observations  on  the  AIO-AFSAT  radio 
link.  Estimates  made  of  the  intensity  scintillation  parameter,  S4,  from  the  SS1ES  observations 
matched  the  observed  S4  derived  from  the  AIO-AFSAT  intensity  data  well,  doing  better  than  the 
WBMOD  climatological  estimates  generated  throughout  the  pass  analyzed.  These  promising 
indications  that  the  data  may  be  useful  in  the  auroral  region  need  to  be  confirmed  by  similar 
analysis  of  data  collected  in  a  repeat  campaign  conducted  in  January  1990. 

There  is  less  certainty  of  the  usefulness  of  these  data  in  the  nighttime  equatorial  region. 
Intense  scintillation  in  this  region  is  dominated  by  depletion-plume  structures,  so  for  the  data  to 
be  of  operational  use  these  structures  must  either  be  directly  observed  by  the  satellite  passing 
through  them  or  their  presence  must  be  inferred  at  some  level  of  confidence  from  observations  of 
the  ambient  plasma  outside  the  plumes.  A  earlier  study*15*  showed  that  the  DMSP  orbit  is  ill- 
suited  for  direct  observation  of  these  plumes  due  to  the  inclination  of  the  orbit  and  the  elongation 
of  the  plumes  along  the  geomagnetic  meridian.  Because  of  this,  there  is  a  significant  probability 
that  the  satellite  may  pass  between  plume  structures  without  observing  them.  The  current  study 
confirmed  this  finding  and  provided  some  support  to  an  assertion  in  the  previous  study  that  it 
may  be  possible  to  infer  the  presence  of  plumes  based  on  a  threshold  effect  in  the  absolute  level 
of  the  electron  density  at  the  satellite  altitude,  This  study  was  too  limited  to  confirm  this 
assertion  or  to  establish  either  the  threshold  level  or  the  degree  of  confidence  one  might  have  in 
predictions  made  of  plume  presence  based  on  this  threshold.  With  several  years  of  SS1ES  data, 
however,  it  may  be  possible  to  address  this  issue  and  resolve  it.  This  is  planned  for  a  follow-on 
effort. 


No  attempt  was  made  in  this  study  to  assess  the  accuracy  of  the  SSIES  scintillation 
estimates  in  locations  and  local  times  other  than  the  auroral  region  in  the  early  evening  hours  and 
the  equatorial  region  at  night.  Assessments  should  be  made  ut  other  latitudes  and  local  times, 
particularly  in  the  polar  cup  and  at  mid-latitudes,  as  it  is  planned  to  use  C^L  values  calculated 
from  the  SSIES  observations  to  provide  near  real-time  updates  to  the  WBMOD  model117*.  Possi¬ 
ble  sources  of  scintillation  data  for  assessments  in  other  ureas  include  routine  and  campaign 
measurements  made  by  GL  ut  the  Thule  geopolc  stution,  measurements  made  by  the  prototype 
TISS  (Translonospheric  Sounding  System)  receiver  in  the  Shetland  Islands,  and  the  extensive 
HiLut/Polar  BEAR  data  sets.  It  is  hoped  to  include  these  assessments  in  the  follow-on  effort. 

In  all,  the  results  of  this  project  were  generally  positive.  Methods  were  developed  to  use 
the  SSIES  data  to  characterize  transionospheric  scintillation  effects,  methods  which  appear  to 
give  reasonable  results  in  the  auroral  region,  While  the  results  ure  not  so  clear  in  the  cquatoriul 
region,  there  may  be  ways  in  which  the  SSIES  data  may  be  used  indirectly  to  infer  the  presence 
of  plume  structures  which  cuuse  the  highest  levels  of  scintillation  encountered  anywhere  in  the 
world.  Further  work  is  required,  but  the  potential  is  there  for  using  these  dutu  to  improve  the 
specification  of  worldwide  scintillation  levels. 
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Appendix.  Topside  Electron  Density  Model 

The  electron  density  profile  model  used  in  this  study  was  developed  specifically  for  use 
with  the  DMSP  SS1ES  sensor  packagel'3!.  The  underlying  structure  of  the  topside  section  of  the 
model  is  that  of  a  two-component  ionosphere  (0+  and  H+)  in  diffusive  equilibrium,  but  it  has 
been  purumetcrized  to  allow  the  model  to  be  fit  to  non-equilibrium  conditions.  Assuming  charge 
neutrality,  the  electron  density  profile  will  be  identical  to  the  height  variation  of  the  ionospheric 
plasma  density,  which  for  this  model  is  given  by 


Ne(h)  »  Np(h)  - 

<x  =  «q  +  (h-400)  [A-2] 

where  Tp  is  the  plasma  temperature  (T^+Te);  Tpo  is  the  plasma  temperature  at  a  reference 
height  (h0);  N0  (0+)  and  N0  (H+)  are  the  number  densities  of  0+  and  H+  at  the  reference 
height;  oq,  cx1(  and  (J  are  profile  adjustment  parameters;  and  g  and  I  are  integral  functions  given 
by 


t  A—  3  ] 
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where  m+  is  the  mass  of  ll+,  and  R.  is  defined  as  the  ratio  of  N(H-t-)  to  N(0+).  Ion  and  electron 
temperatures  are  obtained  from  a  recent  model  based  on  un  analysis  of  data  from  the  AE-C 
satcllitel ,K>.  The  density  ratio,  R,  is  based  on  the  0+  to  H+  transition  height,  h-j«,  defined  as  the 
height  at  which  N(0+)  -  N(H+).  This  parameter  can  be  either  (1)  calculated  from  observations 
of  N(0+),  N(H+),  and  Tp  from  the  SSIES  RPA  sensor,  or  (2)  obtained  from  an  empirical  model 

*  of  l)T  derived  from  published  analyses  of  topside  profiles  from  the  Alouctte  satellites  and  RPA 
datu  from  0G0-6l"l, 

*  The  topside  model  is  fit  to  the  F2  peak  by  means  of  a  parabolic  layer  taken  from  the  Bent 
profile  modcll,l,l  of  the  form 


^  efH*  |  N0(0+)e_16nI  +  N0  (H  +  )  e~x  |  (A-l) 
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Ne(h)  =  NmF2 


t  A-5  ] 


(HFY 

where  NmF2  and  hmF2  are  the  density  and  height  of  the  F2  layer  peak,  and  Yt  is  the  parabolic 
semi- thickness,  The  Yt  parameter  can  be  either  (1)  estimated  in  the  procedure  used  to  fit  the  r 

profile  to  observations,  or  (2)  obtained  from  the  expressions  used  in  the  original  Bent  model. 

Equations  [A-l]  and  [A-51  arc  fit  together  at  the  height  where  the  plasma  scale  height  calculated 
from  the  two  representations  of  Np  (h)  are  equal.  This  must  be  calculated  iteratively,  but  rarely 
requires  more  that  four  or  five  iterations.  This  height  is  also  used  as  the  reference  height,  h0,  for 
the  parameters  in  Equation  (A-l ). 

For  the  sake  of  providing  a  complete  plasma  density  profile,  the  bottomsidc  section  of  the 
Air  Weather  Service  (AWS)  RBTEC  model  is  used  to  describe  the  height  variation  below  the  F2 
peakl20>J*l,  This  model  uses  three  Chapman-function  layers  to  describe  the  three  main  iono¬ 
spheric  layers  (E,  FI,  and  F2).  The  choice  of  models  for  the  bottomsidc  of  the  profile  has  little 
impact  on  the  present  study,  us  the  irregularity  layer  is  assumed  to  start  cither  at  or  just  slightly 
below  the  F2  peak. 
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